Gravitational wave bursts from cosmic strings 



Thibault Damour^-* and Alexander Vilenkin^^ 
Institut des Hautes Etudes Scientifiques, F-91440 Bures- sur-Yvette, France 
Physics Department, Tufts University, Medford, MA 02155, USA 
(26 April 2000) 

Cusps of cosmic strings emit strong beams of high-frequency gravitational waves (GW). As a con- 
sequence of these beams, the stochastic ensemble of gravitational waves generated by a cosmological 
network of oscillating loops is strongly non Gaussian, and includes occasional sharp bursts that stand 
above the rms GW background. These bursts might be detectable by the planned GW detectors 
LIGO/VIRGO and LISA for string tensions as small as Gfj, ~ 10"^"^. The GW bursts discussed here 
might be accompanied by Gamma Ray Bursts. 
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Cosmic strings are linear topological defects that could 
be formed at a symmetry breaking phase transition in the 
early universe. Strings are predicted in a wide class of 
elementary particle models and can give rise to a variety 
of astrophysical phenomena [0. In particular, oscillat- 
ing loops of string can generate a potentially observable 
gravitational wave (GW) background ranging over many 
decades in frequency. The spectrum of this stochastic 
background has been extensively discussed in the liter- 
ature 1^-01, but until now it has been tacitly assumed 
that the GW background is nearly Gaussian. In this 
paper, we show that the GW background from strings 
is strongly non-Gaussian and includes sharp GW bursts 
(GWB) emanating from cosmic string cusps We shall 
estimate the amplitude, frequency spectrum and rate of 
the bursts, and discuss their detectability by the planned 
GW detectors LIGO/VIRGO and LISA. 

We begin with a brief summary of the relevant string 
properties and evolution j^]. A horizon-size volume at 
any cosmic time t contains a few long strings stretching 
across the volume and a large number of small closed 
loops. The typical length and number density of loops 
formed at time t are approximately given by 
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The exact value of the parameter a in (|l|) is not known. 
We shall assume, following that a is determined by 
the gravitational backreaction, so that a ~ TG/i, where 
r ~ 50 is a numerical coefficient, G is Newton's constant, 
and /X is the string tension, i.e. the mass per unit length 
of the string. The coefhcient F enters the total rate of 
energy loss by gravitational radiation d£/dt TG^^. For 
a loop of invariant length / the oscillation period is 
Ti =1/2 and the lifetime is n - l/TG^i ~ t. 

A substantial part of the radiated energy is emitted 
from near-cusp regions where, for a short period of time, 
the string reaches a speed very close to the speed of light 
Q . Cusps tend to be formed a few times during each os- 
cillation period . Let us estimate the (trace-reversed) 
metric perturbation ft.^^ = /i^^ — hh-q^^ emitted near a 



cusp. Let K^v = r-piiysh^v denote the product of the met- 
ric perturbation by the distance away from the string, 
estimated in the local wave-zone of the oscillating string. 
K^i, is given by a Fourier series whose coefficients are pro- 
portional to the Fourier transform of the stress-energy 
tensor of the string: 

yM^(/jA<) ^ y-i I drda^iiX^'X" - X'f'X"') cxp{-ik.X). 

(2) 

Here, X^(r, cr) represents the string worldsheet, 
parametrized by the conformal coordinates t and a. 
[X = OtX, X' = da-X.] In the direction of emis- 
sion n, = (^jk) runs over the discrete set of values 
47rZ~^m(l, n), where m = 1, 2, . . .. Near a cusp (and only 
near a cusp) the Fourier series giving k^^ is dominated 
by large m values, and can be approximated by a con- 
tinuous Fourier integral. The continuous Fourier compo- 
nent (corresponding to an octave of frequency around /) 
«(/) = \ fmf) = \ f\Jdtexp{2nift)K{t) is then given by 



(3) 



In a conformal gauge, the string motion is given by X^^ 
^iXl{a+) + X'^{a_)), where a±=T±a 



and where X!^ 



is a null 4- vector. Further choosing a "time gauge" (r = 
at" = t ) ensures that the time components of these null 
vectors are equal to 1. A cusp is a point on the worldsheet 
where these two null vectors coincide, say AT^^ = = 
(1, ric). One can estimate the waveform (y) emitted near 
the cusp ( localized, say, at ct+ = cr_ = 0) by replacing 
in the integral (^) X'^ by their local Taylor expansions 
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X^{a±) = X!!^ + l^a± + -X^ai + -X^ai 



(4) 



For a given frequency / ;» Tf' , the integral (g) is 
significant only if the angle 9 between the direction of 
emission n and the "3- velocity" of the cusp ric) is smaller 
than about 6',„ = {Ti\f\y^^^. When 9 < 9^, the inte- 
gral can be explicitly evaluated. After a suitable gauge 
transformation one finds [Ol 
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'{f) = -CGfi{27r\f\r'/'Ai^A''_}, 



(5) 



where C = 47r(12)''/3(3r(l/3))-2 and where the hnear 
polarization tensor is the symmetric tensor product of 
Aj. = Xf^/\X±\'^^'^ . The inverse Fourier transform of 
Eq.(p|) gives a time-domain waveform proportional to |t — 
tc\^/ , where tc corresponds to the peak of the burst 
| p2[ . The sharp spike at t = exists only if 6 — 
(i.e. if one observes it exactly in the direction defined 
by the cusp velocity). When 9 1 the spike is 
smoothed over |t — ^ O^Ti. In the Fourier domain 
this smoothing corresponds to an exponential decay for 
frequencies |/| > l/{9^Ti). 

Eq.(^ gives the waveform in the local wave-zone of 
the oscillating loop: h^i, — K^^/rphys- To take into 
account the subsequent propagation of this wave over 
cosmological distances, until it reaches us, one must in- 
troduce three modifications in this waveform: (i) re- 
place J'phys by Qqt where r is the comoving radial co- 
ordinate in a Friedman universe (taken to be flat: ds^ — 
~dt^ + a{tY{dr'^ + r'^dVi^)) and oq = a(to) the present 
scale factor , (ii) express the locally emitted frequency in 
terms of the observed one /cm = (1 + ^)/obs where z is 
the redshift of the source, and (iii) transport the polar- 
ization tensor of the wave by parallel propagation (pp) 
along the null geodesic followed by the GW: 



V(/) = «:^S((l + ^)/)/(«or). 



(6) 



Here, and henceforth, / > denotes the observed fre- 
quency. In terms of the redshift we have a^r = 3io(l — 
(1 -I- z)^^/^), where io is the present age of the universe 
(this relation holds during the matter era, and can be 
used for the present purpose in the radiation era because 
a^jT has a finite limit for large z). 

For our order-of-magnitude estimates we shall as- 
sume that \X±\ ^ 2tt/1. The various numerical factors 
in the equations above nearly compensate each other 
to give the following simple estimate for the observed 
waveform in the frequency domain {h{f) = \f\h{f) = 
|/|/diexp(27ri/t)/,,(t)) 
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Here the explicit redshift dependence is a convenient sim- 
plification of the exact one given above. This result 
holds only if, for a given observed frequency /, the angle 
= cos~"'^(n • ric) satisfies 



(8) 



To know the full dependence of h{f) on the redshift we 
need to express I ^ at m terms of z. We write 

l^atoviiz); ^i{z) = {l + z)-^/\l + z/z,q)-^/^. (9) 

Here Zeq — 2.4 x lO^'iloft-o — lO'^-^ is the redshift of equal 
matter and radiation densities, and we found it conve- 
nient to define the function fi{z) which interpolates be- 
tween the different functional z-dependences of I in the 



matter era, and the radiation era. [We shall systemat- 
ically introduce such interpolating functions of z, valid 
for all redshifts, in the following.] Inserting Eq.(||) into 
Eq.(|) yields 

h{f,z)^Gfia^/^{fto)-^/^Mz), 
^h{z) = z-\\ + z)-i/3(l + z/z,,)-i/3. (10) 

We can estimate the rate of GWBs originating at cusps 
in the redshift interval dz, and observed around the fre- 
quency /, as dN ~ \0'^{l + z)~'^v{z)dV{z). Here, the 
first factor is the beaming fraction within the cone of 
maximal angle 9m{f,z), Eq.(||), the second factor comes 
from the relation diobs = (1 + z)dt, v{t) ~ fcn'i{t)/Ti 
2a~^t~^ is the number of cusp events per unit spacetime 
volume, /c ~ 1 is the average number of cusps per oscilla- 
tion period of a loop, T/ ^ at/2, and dV{z) is the proper 
volume between redshifts z and z-l-dz. In the matter era 
dV = 547ri3[(l -I- z)i/2 _ i]2(i + ^)-ii/2^2^ ^i^jig tj^g 

radiation era dV = 727rij](l -I- Ze^)i/2(1 + z)-^dz. The 
function N[J, z) = dN/dhiz can be approximately rep- 
resented by the following interpolating function of z 



7V(/,z) ^ 10%-ia-«/3(/to)-'/Vn(^), 
^„(z) = z3(l-Hz)-^/6(l-t-z/ze,)"/6. 



(11) 



The observationally most relevant question is: what 
is the typical amplitude of bursts /i^"'''*(/)that we can 

expect to detect at some given rate N, say, one per year? 
Using TV = /q" 7V(/, z)dlnz ^ N{f,Zm), where Zm is 
the largest redshift contributin g to TV, one can estimate 



^burst^j-^ by solving for z in Eq.(|llD and substituting the 



N 



result z = Zm{N,f) in Eq.(p^). The final answer has a 
different functional form depending on the magnitude of 
the quantity 



y{NJ)^ 10-' Ntoa""^' if tof/\ 



(12) 



Indeed, if ?/ < 1 the dominant redshift will be Zm{y) < 1; 
while, ii 1 < y < Zeq^^^^, 1 < z„i(y) < Zg^, and if y > 
Zeq'^^^^, Zm{y) > Zeq- We cau again introduce a suitable 
interpolating function g{y) to represent the final result 
as an explicit function of y: 



^ burst 
''N 



{f)^Gtia'/\ft,r^/'g[y{NJ)], 



9{y) = y-i/3(i + 2;)-i3/33(i + j,/(,^^)n/6)3/ii. (13) 

The prediction Eq. ( p^ ) for the amplitude of the GWBs 
generated at cusps of cosmic strings is the main new re- 
sult of this work. To see whether or not these bursts can 
be distinguished from the stochastic gravitational wave 
background, we have to compare the burst amplitude 
(|l3| ) to the rms amplitude of the background, /irms, at 
the same frequency. We define h^ms as the "confusion" 
part of the ensemble of bursts Eq.(|l3|), i.e. the super- 
position of all the "overlapping" bursts, those whose oc- 
currence rate is higher than their typical frequency. This 
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is estimated by multiplying the square of Eq. (|T^) by the 
number of overlapping bursts within a frequency octave 
Uzif) = f^^N{f,z) , and then integrating over all Inz 
such that nz{f) > 1 , and Om{f, z) < 1: 

hLsif) = J h\j,z)nz{J)d\nzH{nz-l)H{l-e^), 

(14) 

where H denotes Heaviside's step function. Eq.(jlJ) dif- 
fers from previous estimates of the stochastic background 
(beyond the fact that we use the simplified loop 
density model Eq.(|l])) in that the latter did not incor- 
porate the restriction to nz{f) > 1, i.e. they included 
non-overlapping bursts in the average of the squared GW 
amplitude. 

It is easily checked from Eq.(|l|) that /i^""''* IS a mono- 
tonically decreasing function of both N and /. These de- 
cays can be described by (approximate) power laws , with 
an index which depends on the relevant range of domi- 
nant redshifts; e.g., as N increases, h^^'^^*^ decreases first 
like iV"^/^ ( in the range < 1), then like iV~^/^^ (when 
1 < Zm < Zeq), and finally like N^^^^^ (when Zm > Zeq). 
For the frequency dependence of h^^^'^^ , the correspond- 
ing power-law indices are successively: —5/9, —9/11 and 
—7/11. [These slopes come from combining the basic 
j-i/3 dependence of the spectrum of each burst with 
the indirect dependence on / of the dominant redshift 
Zmia^N, /).] By contrast, when using our assumed link 
G/i ~ a/50 between the string tension fj, and the pa- 
rameter a, one finds that the index of the power-law de- 
pendence of /jburst ^pQ]-^ Q, takes successively the values 
-1-7/9,-3/11 and -1-5/11. Therefore, in a certain range 
of values of a (corresponding to 1 < Zm{a,N, f) < z^q) 
the GWB amplitude (paradoxically) increases as one de- 
creases a , i.e. G/i. 

In Fig. 1 we plot (as a solid line) the logarithm of 
the GW burst amplitude, log;^o(''''^^'^'^*)' ^ function of 
logiQ{a), for TV = 1 yr~^, and for f — f^. — 150 Hz. 
This central frequency is the optimal one for the detec- 
tion of a /"^/'^-spectrum burst by LIGO. We indicate 
on the same plot (as horizontal dashed lines) the (one 
sigma) noise levels of LIGO 1 (the initial detec- 

tor), and LIGO 2 (its planned advanced configuration). 
The VIRGO detector has essentially the same noise level 
as LIGO 1 for the GW bursts considered here. These 
noise levels are defined so that the integrated optimal 
( with a matched filter oc |/|~^^'^) signal to noise ra- 
tio (SNR) for each detector is SNR = /i'^"'^"*(/c)//i''°'"^ 
The short-dashed line in the lower right corner is the 
rms GW amplitude, Eq. ([T^) . One sees that the burst 
amp litude stands well above the stochastic background 
0. Clearly the search by LIGO/VIRGO of the type of 
GW bursts discussed here is a sensitive probe of the ex- 
istence of cosmic strings in a larger range of values of a 
than the usually considered search for a stochastic GW 
background. 
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FIG. 1. Gravitational wave amplitude of bursts emitted by 
cosmic string cusps in the LIGO/VIRGO frequency band, as 
a function of the parameter a = 50G/i. [In a base-10 log-log 
plot.] The horizontal dashed lines indicate the one sigma noise 
levels (after optimal filtering) of LIGO 1 (initial detector) and 
LIGO 2 (advanced configuration) .The short-dashed line indi- 
cates the rms amplitude of the stochastic GW background. 

From Fig. 1 we see that the discovery potential of 
ground-based GW interferometric detectors is richer than 
hitherto envisaged, as it could detect cosmic strings in 
the range a ^ 10~^°, i.e. G/z ^ 10~^^. Let us also note 
that the value of a suggested by the (superconducting-) 
cosmic-string Gamma Ray Burst (GRB) model of ReL 
[]l3| , namely a ~ 10~®, nearly corresponds, in Fig. 1, 
to a local maximum of the GW burst amplitude. [This 
local maximum corresponds to ~ 1- The local min- 
imum on its right corresponds to Zm ~ Zeq.] In view 
of the crudencss of our estimates, it is quite possible 
that LIGO 1 /VIRGO might be sensitive enough to detect 
these GW bursts. Indeed, if one searches for GW bursts 
which are (nearly) coincident with (some [^) GRB the 
needed threshold for a convincing coincident detection is 
much closer to unity than in a blind search. [ In a blind 
search, by two detectors, one probably needs SNRs ~ 4.4 
to allow for the many possible arrival times. Note that 
the optimal filter, /jtcmpiatc^j) ^ |j|-i/3^ fo^, ^^j. gwBs 

is paramctcr-free.] 

In Fig. 2 we plot log]^Q(/i'^"''''*) as a function of log]^Q(a) 
for iV = 1 yr-i, and for / = = 3.9 x IQ-^ Hz. 
This frequency is the optimal one for the detection of 
a /^^/'^ GWB by the planned spaceborne GW detec- 
tor LISA. [In determining the optimal SNR in LISA we 
combined the latest estimate of the instrumental noise 
p6[ with estimates of the galactic confusion noise [p7[|.] 
Fig. 2 compares /i''"''"*(/c) to both EISA's (fifiered) noise 
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level and to the cosmic-string-generated stochastic 

background /i""'', Eq.(14). The main differences with the 
previous plot are: (i) the signal strength, and the SNR, 
are typically much higher for LISA than for LIGO, (ii) 
though the GW burst signal still stands out well above 
the rms background, the latter is now higher than the 
(broad-band) detector noise in a wide range of values 
of a. LISA is clearly a very sensitive probe of cosmic 
strings. It might detect GWBs for values of a as small 
as ~ 10~^^'^. [Again, a search in coincidence with GRBs 
would ease detection. Note, however, that, thanks to 
the lower frequency range, even a blind search by the 
(roughly) two independent arms of LISA would need a 
lower threshold, ^ 3, than LIGO.] 



space density of loops. This estimate may be accurate in 
the matter era but is probably too small by a factor ~ 10 
in the radiation era . Taking into account this increase 
of ni would reinforce our conclusions in making easier (in 
some parameter range) the detection of GWBs.(iv) Fi- 
nally, we disregarded the possibility of a nonzero cosmo- 
logical constant which would introduce some quantitative 
changes in our estimates. 

The work of A.V. was supported in part by the Na- 
tional Science Foundation. 



FIG. 2. Gravitational wave amplitude of bursts emitted by 
cosmic string cusps in the LISA frequency band, as a func- 
tion of the parameter a — bOGfi. [In a base- 10 log- log plot.] 
The short-dashed curve indicates the rms amplitude of the 
stochastic GW background. The lower long-dashed line in- 
dicates the one sigma noise level (after optimal filtering) of 
LISA. 

We shall discuss elsewhere the consequences for the 
interpretation of the pulsar timing experiments of the 
GWB-induced non-Gaussianity of the stochastic GW 
background [0. 

When comparing our results with observations, one 
should keep in mind that the model we used for cosmic 
strings involves a number of simplifying assumptions, (i) 
All loops at time t were assumed to have length / ~ at 
with a ~ FG/Lt. It is possible, however, that the loops 
have a broad length distribution n{l,t) and that the pa- 
rameter a characterizing the typical loop length is in the 
range TGji < a ^ 10^"^. (Here, the upper bound comes 
from numerical simulations and the lower bound is due 
to the gravitational backreaction.) (ii) We also assumed 
that the loops are characterized by a single length scale 
I, with no wiggliness on smaller scales. Short-wavelength 
wiggles on scales <C TGfjd are damped by gravitational 
backreaction, but some residual wiggliness may survive. 
As a result, the amplitude and the angular distribution of 
gravitational radiation from cusps may be modified. (Hi) 
We assumed the simple, uniform estimate Eq.(nh for the 
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